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Two recent theoretical works studied the role of Kitaev interactions in the newly observed incommensurate 
magnetic order in the hyper-honeycomb (J-LUIrCE) and stripy-honeycomb ( 7 -Li 2 lr 03 ) iridates. Each of these 
works analyzed a different model (JAT versus coupled zigzag chain model) using a contrasting method (classi¬ 
cal versus soft-spin analysis). The lack of commonality between these works precludes meaningful comparisons 
and a proper understanding of these unusual orderings. In this study, we complete the unfinished picture initiated 
by these two works by solving both models with both approaches for both 3D honeycomb iridates. Through 
comparisons between all combinations of models, techniques, and materials, we find that the bond-isotropic 
JAT model consistently predicts the experimental phase of /ALUMT regardless of the method used, while the 
experimental phase of y-LEIrCE can be generated by the soft-spin approach with eigenmode mixing irrespec¬ 
tive of the model used. To gain further insights, we solve a ID quantum spin-chain model related to both 3D 
models using the density matrix renormalization group method to form a benchmark. We discover that in the 
ID model, incommensurate correlations in the classical and soft-spin analysis survive in the quantum limit only 
in the presence of the symmetric-off-diagonal exchange T found in the JAT model. The relevance of these 
results to the real materials are also discussed. 


I. INTRODUCTION 

Kitaev’s honeycomb model, which hosts an exact Z2 spin 
liquid 1 , has been identified to play a crucial role in the 
low-energy description of the quasi-two-dimensional lay¬ 
ered honeycomb iridates A 2 Ir03 (/l=Li, Na). 2 : A flurry of 
theoretical 5 " 12 and experimental 1 ' 20 studies have since ex¬ 
amined the material’s properties. Most recently, two three- 
dimensional ( 3 D) analogues of these layered materials have 
been synthesized 21-22 , giving hope to being the first realiza¬ 
tion of a 3 D spin liquid. These discoveries have spurred many 
theoretical studies on the nature of these spin-orbit coupled 
Mott insulators. 23 " 31 However, much like their 2 D analogues, 
these two materials—/3-Li 2 Ir03 32 and y-LEh'O.a ' (hereafter 
/ 3 -LIO and 7-LIO)—are magnetically ordered rather than ex¬ 
hibiting spin-liquid behavior. To better understand the role 
of Kitaev physics in these systems, their magnetic behav¬ 
ior must first be scrutinized. As of late, two theoretical 
studies have attacked the problem from two complementary 
directions: 23-29 each of these studies proposed a pseudospin 
1/2 model describing the low-energy physics of the two iri¬ 
dates and analyzed their respective model using a different 
approach. The first work examined the JAT model using a 
classical approach 28 , while the second analyzed the coupled 
zigzag chain (CZC) model using a soft-spin method 29 (see 
Sec. Ill for definition of the models). The goal of both works 
is to utilize the detailed experimental magnetic structure to 
constrain their respective model to the region of the parameter 
space in which J-LUIrO.j and 7-Li 2 Ir03 reside. One common 
conclusion between the two studies is that the dominance of 
a ferromagnetic Kitaev term and the presence of subdominant 
exchange terms (including the antiferromagnetic Heisenberg 
exchange) generate spiral orders through frustration arising 


from anisotropic exchanges. 

While both works unearthed certain features of the mag¬ 
netism in these new iridates, their approaches and results have 
important differences. Without a reliable method of solving a 
frustrated quantum spin model in three dimensions, a defini¬ 
tive evaluation of these two models and methods is beyond 
reach. All hope is not lost, however: a careful comparison of 
the models and techniques used can offer a much more co¬ 
hesive and encompassing understanding of the magnetism in 
these 3 D honeycomb iridates. 

In this study, we complete the picture depicted in these pre¬ 
vious works by examining the JKT and CZC models using 
both classical and soft-spin approaches. This work compares 
all combinations of models and methods, culminating in an 
understanding greater than the sum of its parts. Our main re¬ 
sult is summarized in Table I, where we compare the resulting 
ground states of all eight combinations of materials, methods, 
and models. With these results, we conclude that the JKT 
model on the hyper-honeycomb is most robust as it predicts 
the experimental phase under all methods. Additionally, we 
find that the 7-LIO magnetic structure can only be captured 
with the soft-spin approach with linear combination of multi¬ 
ple eigenmodes in both the JKT and CZC models. 

Ultimately, to understand how each of these methods suc¬ 
ceeds or fails to capture the quantum ground state properties, 
a comparison with a quantum treatment would be highly de¬ 
sirable. However, these frustrated quantum mechanical spin 
models are difficult to solve in three dimensions. Hence, we 
analyze a ID spin-chain model that is related to both the JKT 
and CZC models. In one dimension, the quantum solution 
of the spin model was obtained numerically using the den¬ 
sity matrix renormalization group (DMRG) technique. We 
compare the classical solution and soft-spin analysis with the 
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Lattice 

Experimental 

Model 

Classical 

Soft-spin 




no linear combination 

with linear combination 

/3-LIO 

(A a , C b , F c ) 

JKT 

(A a ,C b ,F c ) 

(A a , C b , F c ) 

(A a , C b , F c )“ 

CZC 

Commensurate 

(A a ,0 ,F C ) 

(A a ,C b ,F c ) 

7-LIO 

(A a ,F b ,F c ) 

JKT 

(A a , Cb, F c ) 

(Aa, Cb, F c ) 

(A a ,F b ,F c ) 

CZC 

Commensurate 

(A a ,0,F c ) 

(A a ,F b ,F c ) 


a Since the symmetry of the experimental spiral phase was already obtained without linear combination of multiple eigenmodes, no such mixing was 
performed. 


TABLE I. Summary of results for all combinations of model, lattice, and method. The irreducible representations of the a, b, and c components 
of the moments in the ground states are given in a shorthand notation. The entry is in bold if the irreducible representations correspond to the 
experimental phase, i.e. the theoretical phase matches in symmetry with the observed ordering. 




FIG. 1. (Color online) Lattices of the hyper- and stripy-honeycomb. 
The red and green bonds are the symmetry-equivalent X and Y 
bonds. The blue bonds are the Z bonds, which are symmetry- 
inequivalent to the X/Y bonds. For the stripy-honeycomb lattice, 
there are two inquivalent sets of Z bonds. The bonds between black 
sites are the Z\ bonds while the bonds between red sites are the Z 2 
bonds. The labeling of these bonds correspond to the Ising compo¬ 
nent present in the definition of the Kitaev interactions in both the 
JAT and CZC models. The numbers indicate the sublattice labeling 
used in this work. 


quantum results and find that incommensurate correlations 


found in both classical and soft-spin methods only persist in 
the quantum limit when the symmetric off-diagonal exchange 
r is finite. Although useful as an illustrative tool, we also cau¬ 
tion the extrapolation of these results to higher dimensions, 
noting that quantum models can behave drastically differently 
as the dimensionality of the system changes. 

In Sec. II, we briefly provide background information on 
the two iridates in question: the hyper-honeycomb /3-LIO and 
the stripy-honeycomb 7 -LIO. In Sec. Ill, we outline the differ¬ 
ences between the two spin models—the JKT and the CZC 
models. Then, in Sec. IV, we delve into the details of the 
classical and soft-spin approaches, applying these methods to 
obtain phase diagrams for the two models at hand. We con¬ 
clude this section with a discussion about the implications of 
our findings. We present our ID spin chain results in Sec. V, 
and relate these results to those obtained earlier for the 3D 
models. Lastly in Sec. VI, we provide a summary of our re¬ 
sults and an outlook on the topic of magnetism in these 3D 
iridate materials. 


II. TWO IRIDATES 

The two 3D honeycomb iridates are structural variants 
(polymorphs) of Li 2 IrC> 3 , and the elementary building blocks 
of both compounds is the IrOr; octahedron. Each octahedron 
shares three edges with its three neighbors, defining a tri- 
coordinated network of Ir ions. The resulting networks differ 
between the two compounds, as seen in Fig. 1. Distinguish¬ 
ing features include the number of symmetry-inequivalent 
nearest-neighbor (NN) bonds and the symmetries they pos¬ 
sess. In the hyper-honeycomb, there are only two symmetry- 
inequivalent sets of bonds: the X/Y bonds that possesses in¬ 
version symmetry and the Z bonds that possesses 222 point 
group symmetry. In contrast, there are three symmetry- 
inequivalent sets of NN bonds in the stripy-honeycomb: the 
X/Y bonds that do not possess any symmetry, the Z\ bonds 
that possess 2/m symmetry, and the Z 2 bonds that possess 
222 symmetry. These distinctions are important in the con¬ 
struction of models for these materials and in results discussed 
in Sec. IV C . 34 

The magnetism found in these two iridates have many strik¬ 
ing similarities. Both materials undergo a magnetic phase 
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transition at T ~ 38 I \. 21,22 The magnetic order in both cases 
is a non-coplanar, incommensurate spiral order, where mo¬ 
ments on neighboring sublattices rotate in the opposite sense 
(also known as counter-rotation of spirals). Moreover, the 
wavevector of both spirals are equal within experimental res¬ 
olution: ~(0.5 7,0,0) in the hkl notation of each compound’s 
respective orthorhombic unit cell/ 2 ' ” 

While very similar, differences in the detailed magnetic 
structure have been identified. ' 2 ' ” In the hyper-honeycomb 
case, the spiral order transforms under a single irreducible 
representation (irrep) of the magnetic space group: T 4 (see 
Appendix A for details )/ 2 In other words, when applied to 
the spiral order, all symmetries of the lattice that leaves the 
ordering wavevector ~ (0.57,0,0) invariant results in a fac¬ 
tor of ± 1 , where the sign is dictated by the character of the 
irreducible representation. In contrast, the spiral in the stripy- 
honeycomb case transforms in a more complicated fashion: 
the moments parallel to the orthorhombic b direction trans¬ 
form under a different irrep (Ty) than the moments perpendic¬ 
ular to the b direction (I^). 3 -’ Alternatively, in terms of mag¬ 
netic basis vectors, this implies that the magnetic structure of 
/3-LIO transforms as (A a . C/,, F c ) while the magnetic struc¬ 
ture of 7 -LIO transforms as (A a , J/, F c ). This subtle differ¬ 
ence between the two spiral orders is outlined in the first two 
columns of Table I and plays a crucial role in our discussion 
when comparing the two models and the two methods. 


III. TWO MODELS 


A low-energy description of these Mott insulators with 
strong SOC can be obtained by starting in the atomic pic¬ 
ture. In this limit, the partially-filled Ir 4+ ions are in the 
5 d 5 configuration. In the presence of large octahedral crystal 
fields and strong SOC, an effective low-energy description can 
be derived from a single electron occupying the high-energy 
j e ff = 1/2 doublet. In the Mott insulating limit, the low- 
energy degrees of freedom of these half-filled doublets can 
be represented by j e ff = 1/2 pseudospins. Due to the strong 
SOC, these pseudospins are expected to interact via strongly 
anisotropic exchanges ( i.e . strong exchange-anisotropy). This 
is supported by perturbative calculations that capture both the 
effects of virtual states in the presence of Hund’s coupling 
and oxygen-mediated superexchange mechanisms. We now 
explore two pseudospin models motivated by these symmetry 
considerations. 


A. JKT model 


The JKV model proposed in Ref. 28 assumes that all 
bonds have the same exchange parameters (bond-isotropy). In 
addition, the model assumes all bonds take on three exchange 
parameters J, K, and T, resulting in the Hamiltonian 


H = 


{ij)£ai 3 ( 7 ) 


JSi ■ s, + Ks?sn + e^rs?s: 


(i) 


where the pseudospin at site i is denoted as S, , NN bonds 
are denoted as (ij), and each NN bond is labelled by 7 € 
(X, Y,Z). The shorthand (ij) £ a/. 3 ( 7 ) means (ij) £ 7 and 
an implicit sum is taken over only a and f3 where a//? / 7 . 
The phase angle 0 r) £ {0, 7 r} mod 27 t is bond-dependent and 
is partially constrained by the symmetry of the lattice. In the 
model considered in Ref. 28 and here, it is given by 6ij = 
Tt(fij ■ v + 1) where fij is the unit vector from site i to j, 

u=^|(| + |), and b, c are the conventional lattice vectors. 

B. Coupled zigzag chain model 

The coupled zigzag chain (CZC) model proposed in Ref. 
29 assumes that the X/Y bonds and Z are distinct (for 7 -LIO, 
Z -2 and Z\ bonds are assumed to possess the same exchange 
parameters). Hence, this is an intrinsically bond-anisotropic 
model. Furthermore, the number of exchange parameters 
have been restricted: while the Z bonds (Z 2 and Z\ bonds 
in 7 -LIO) are assumed to have three parameters similar to the 
JKT model, the X/Y bonds possess only two exchanges that 
parametrize the Heisenberg and Kitaev exchanges. These ex¬ 
changes are related as the X/Y bonds and Z bonds share the 
same J and K values via 

H= J2 [JS i -S j +KS/S]]+ J2 Ur i:j -SiXfij-Sj). 
(ij) 67 (ij)£Z 

In the above, (ij) denotes NN sites, ./. K. I c are the Heisen¬ 
berg, Kitaev, and Ising couplings respectively, 7 is the Kitaev 
component and label on bond ij, r 1/} is the unit vector from 
site i to j, and the Ising term only sums over the Z bonds. 
One can also re-write the exchanges in terms of the JKT 
parametrization on each bond to manifestly show the bond- 
anisotropy 

Jz = J + \h I<z = K- h c T z = l -I c 

Jx/Y = J K x /y = K r x/Y = 0, 

where the bond label of each exchange is indicated by the sub¬ 
script. 


IV. TWO APPROACHES 


In the two studies Ref. 28 and 29, not only were the mod¬ 
els considered different, the approaches used to identify the 
ground states of each model also differed. Here we first pro¬ 
vide details on both methods then subsequently complete the 
analysis of the prior studies by applying both methods to both 
models on the two iridates in consideration. 


A. Classical analysis 

In the classical limit, the quantum mechanical pseudospins 
are treated as constant length vectors. For S = 1/2, the clas¬ 
sical treatment is equivalent to a variational method where the 
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(a) Hyper-honeycomb: classical 
analysis 


(b) Hyper-honeycomb: soft-spin 
analysis 


(a) Hyper- and stripy-honeycomb: (b) Hyper- and stripy-honeycomb: 
classical analysis soft-spin analysis 



(c) Stripy-honeycomb: classical 
analysis 


id) Stripy-honeycomb: soft-spin 
analysis 


FIG. 2. (Color online) JKT model: comparing the classical 
and soft-spin (with no mixing) approaches on the hyper-honeycomb 
and stripy-honeycomb lattices. The classical approach stipulates the 
constant-spin-length constraint whereas the soft-spin approach does 
not. LT fails in the region bounded by the red dotted line. The non- 
coplanar spiral (NCsp a ) state exist in all phase diagrams—it is the 
experimentally observed phase in the hyper-honeycomb /3-LIO and 
closely related to the experimental phase in the stripy-honeycomb 
7 -LIO. The values at the circle boundary are 0-values, whereas val¬ 
ues inside are ^-values; see Eq. 2 and surrounding main text for 
parametrization used. 


ansatz is restricted to a site-factorized product state. This is 
because the product state ansatz must satisfy the classical unit 
length constraint due to the unique Bloch sphere nature of 
S' = 1/2 spins: the expectation value (S) for any pure state 
forms a constant length vector. As such, to capture long-range 
ordered states, the classical approach is a useful first method 
to employ. 

To solve the classical problem (or equivalently, the varia¬ 
tional problem), we have various tools at our disposal. The 
Luttinger-Tisza (LT) approximation is an important method 
that serves to provide a lower bound to the classical ground 
state energy and it also identifies the exact classical ground 
state when the method succeeds. Success is debited as find- 
ing a solution that satisbes the constant spin-length constraint 
across all sites and is typically met in regions of parameter 
space where geometrical and/or exchange frustration is mini¬ 
mal. 


FIG. 3. (Color online) CZC model: comparing the classical and 
soft-spin (with no mixing) approaches on the hyper-honeycomb and 
stripy-honeycomb lattices. The hyper- and stripy-honeycomb results 
agree within each method hence only one figure is shown per method. 
LT fails in the region bounded by the red dotted line. The classical 
method finds commensurate ordering in the whole parameter space 
considered, whereas the soft-spin approach finds a coplanar spiral in 
the LT-failed region. See text for additional details on phases present 
and discussion on the comparison. 


Procedurally, the essence of the LT method involves tak¬ 
ing a Fourier transform of the real-space exchange Hamil¬ 
tonian into momentum-space and diagonalizing the resultant 
matrix at all momenta to yield an LT band structure. Then, we 
identify the momenta that globally minimize such band struc¬ 
ture and examine the corresponding degenerate eigenspace. 
Lastly, we attempt to construct an eigenvector within the de¬ 
generate eigenspace that, when Fourier transformed back into 
real-space, yields a constant spin-length configuration. If such 
a configuration is found, then the state is a classical ground 
state and the energy corresponds to the state’s eigenvalue. If 
no such configuration can be found, the LT method is deemed 
to have failed and the lowest eigenvalue is the lower bound to 
the classical ground state energy. In the regions of parameter 
space where LT fails, we in turn employ numerical methods 
such as simulated annealing to find the classical ground state, 
supplemented with the knowledge of the lower bound in the 
classical energy from the LT analysis. 


1. Previous results from JKT model 

We first review the classical results of the JKT model. In 
Fig. 2a and 2c, we reproduce the classical phase diagram of 
the JKT model for both 3D honeycomb systems from Ref. 
28 with J > 0, K < 0, and I’ < 0, which contains the experi¬ 
mentally relevant spiral phases. The exchange interactions are 
parametrized in the polar plot as 

(J,iv,r) = (sin 6 cos 0, sin 8 sin 0, cos 6), (2) 

where 0 is the angular coordinate and r = (7r — 0) is the ra¬ 
dial coordinate. As such, the outer edge of the quarter-circle 
is the Heisenberg-Kitaev limit, the left edge is the J = 0 limit. 
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and the top edge is the K = 0 limit. Outside the red dotted 
lines, LT succeeds and an antiferromagnetic phase (AF a ) ex¬ 
ists. Within the red dotted lines, LT fails and we subsequently 
employed simulated annealing to identify the ground states in 
this outlined region. 

Hyper-honeycomb: The NCsp a (A'on-Coplanar spiral) 
phase in the dotted region reproduces the experimentally ob¬ 
served magnetic order. It is a non-coplanar, counter-rotating 
spiral order with the experimentally observed broken sym¬ 
metries and an ordering wavevector in the hOO direction. In 
terms of the magnetic basis vectors, the NCsp a phase trans¬ 
forms as (A a ,Cb, F c ), i.e. as a single-irrep T 4 . The order¬ 
ing wavevector in the NCsp a region continuously changes as 
a function of J, K , and T, and the experimental wavevector 
is contained within. In addition to the NCsp a phase, there is 
also another non-coplanar spiral with wavevector in the OfcO 
direction (NCsp6) and a skew-stripy phase (SS/,). 

Stripy-honeycomb: In contrast to the hyper-honeycomb 
case, none of the phases in the stripy-honeycomb phase dia¬ 
gram reproduces the experimental phase exactly, although the 
NCsp a phase in the stripy-honeycomb model is similar to the 
experimental phase. In particular, the NCsp Q phase transforms 
as the single irrep 1?4 with symmetry {A a , Cb, F c ), which dif¬ 
fers from the experimental phase in the b component. Like the 
hyper-honeycomb case, the wavevectors in the NCsp a region 
continuously change as a function of the exchange parameters 
and the region contains the experimental wavevector. 

For further details on the classical results for the JKT 
model including the existence of other magnetic orders and 
other parameter regimes, we refer the reader to Ref. 28. 


2. New results from CZC model 

In Fig. 3a, we applied the classical analysis on the CZC 
model for both the hyper-honeycomb and stripy-honeycomb 
lattices in the same parameter region as that of Ref. 29, i.e. 
0 < J < 0.5, K = — 1, and 0 > I c > —1.0. The phase di¬ 
agrams of both the hyper- and stripy-honeycomb CZC model 
yield identical phase boundaries and hence only one phase di¬ 
agram is shown. Outside the red dotted region where either 
the Heisenberg exchange (J) or the Ising exchange {If) dom¬ 
inates, LT succeeds and identifies two distinct commensurate 
phases. In the large J/|A'| region, the ground state is a skew- 
stripy phase (SS/J, while in the large IJK region, a collinear 
ferromagnetic ground state (FM C ) in the c direction is found. 

In the dotted region where LT fails, simulated annealing 
identifies commensurate ground states. The FM C phase that 
was found in the large I c /K region enlarges into the dotted 
region, and a new skew-stripy phase (SS a /b) emerges between 
the FM C and SSf, phase. The experimental spiral phase does 
not appear in the classical analysis of the CZC model in either 
crystal structures. 


B. Soft-spin analysis 

In the soft-spin analysis as proposed in Ref. 29, one initially 
follows the same procedure as the LT method as outlined in 
Sec. IV A, but disregards the constant spin-length constraint. 
It was suggested that the minimum eigenvalue solution, which 
may violate the spin-length constraint, can be considered as a 
soft-spin solution. 29 If the constant spin-length condition is 
satisfied, however, the solution will be the same as the clas¬ 
sical solution. It was further proposed that one can construct 
a more complex spin structure with the same wavevector by 
forming a linear combination, or mixing, of the lowest and cer¬ 
tain higher-eigenvalue eigenvectors, which may be considered 
as a potential candidate solution for the quantum model. 24 
Here we will follow their prescription and investigate the re¬ 
sults in both models. 


1. Previous results from CZC model 

In Fig. 3b, we reproduced the soft-spin phase diagram ob¬ 
tained in Ref. 29, which applies to both hyper- and stripy- 
honeycomb CZC models. As explained in Sec. IV B, the dif¬ 
ference between the classical and soft-spin results lies within 
the dotted region where LT fails. In this region, the soft-spin 
analysis finds the Csp Q phase with a continuously changing 
wavevector as a function of J/\K\ and I c /K. This phase 
is a coplanar, counter-rotating spiral with vanishing moments 
along the orthorhombic b direction and a wavevector along the 
IiOO direction. 

This state differs from the experimental phase because mo¬ 
ments are only in the a-c plane, hence yielding the symme¬ 
try (A a ,0, F c ) (0 is used to denote the vanishing b compo¬ 
nent). To obtain a finite component along the b direction 
while ensuring the symmetry of the experimental phase, the 
eigenvector of the fifth-lowest eigenvalue was linearly com¬ 
bined with the eigenvectors of the lowest eigenvalue. 24 This 
mixed state has the symmetry (A a , Cb, F c ) in the case of the 
hyper-honeycomb, and ( A a ,Fb, F c ) in the case of the stripy- 
honeycomb, in agreement with experimental results. 29 It was 
proposed that this mixture of eigenvectors would result in a 
candidate solution for the quantum model. 24 


2. New results from JKF model 

In Fig. 2b and 2d, we have computed the phase diagram 
of the JKT model using the soft-spin approach. For both 
the hyper-honeycomb and stripy-honeycomb cases, in the re¬ 
gion where one would find the NCsp a phase using the clas¬ 
sical approach, we find phases with the same broken sym¬ 
metry albeit with modified ordering wavevector lengths and 
non-constant spin-lengths. In this region of both lattices, the 
ordering wavevector continuously changes and contains the 
experimentally verified wavevector. 

Hyper-honeycomb: The symmetry of this lowest- 
eigenvalue state is (A a ,Cb, F c ), and since it agrees with 
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the experimental phase, a linear combination with higher- 
eigenvalue eigenvectors is not necessary. In other words, the 
JKT model on the hyper-honeycomnb lattice succeeds in 
finding the experimental phase in both the classical analysis 
and the soft-spin analysis (without mixture). In comparison 
with the classical results, both the NCsp a and NCsp/, phases 
enlarge and meets the boundary of where LT fails. The SSf, 
phase that was found in the classical analysis does not appear 
in the soft-spin analysis. 

Stripy-honeycomb: The lowest-eigenvalue state has sym¬ 

metry (A ai C'b , F c ), which is identical to the NCsp a phase in 
the classical analysis. A state with the experimental symme¬ 
try {A a , Ft, F c ) can be obtained if higher-eigenvalue eigen¬ 
vectors are mixed, similar to the results outlined above for 
the CZC model. In addition to this phase, there are also 
two additional non-coplanar spirals, NCsp', and NCsp/,, and a 
non-coplanar, commensurate, antiferromagnetic phase, AF ac . 
NCsp', has an ordering wavevector along the /;00 direction 
but has a different symmetry from NCsp a . On the other hand, 
NCsp/, has an ordering wavevector along the OfcO direction in 
analogy to the similarly named phase in the hyper-honeycomb 
phase diagram. The non-coplanar, commensurate antiferro¬ 
magnetic phase AF ac lies close to the Kitaev limit and is dis¬ 
tinct from the SS x / y phase predicted in the classical analysis. 
This phase consists of moments aligned antiferromagnetically 
in the c direction on the Z\ bonds, and antiferromagnetically 
in the a direction on the Z 2 bonds with a small ferromagnetic 
component along the b axis. Overall, the phase does not have a 
net moment as the ferromagnetic tilts on the various Z 2 bonds 
within the unit cell are anti-aligned. 


C. Discussion 

Table I summarizes our main results and allows for an easy 
comparison between all methods, models, and lattices at a 
glance. The table lists the symmetry of the experimentally 
observed ordering for each lattice, followed by the symmetry 
of the generated phases present for each model using classical 
analysis and the soft-spin analysis with and without mixing of 
eigenmodes. 

In the /3-LIO case, the JKT model generates the experi¬ 
mental phase under all methods in consideration. In contrast 
to the robustness of the JKT model, the CZC model can only 
generate the experimental phase under the soft-spin approach 
with mixing of eigenmodes. Since the JKT model is mani¬ 
festly bond-isotropic while the CZC model is not, this differ¬ 
ence in robustness in generating the experimental phase may 
be attributed to the bond-isotropy of the real material. The 
bond-isotropy of /3-LIO is supported by the near-ideal Ir-Ir 
and Ir-O bond lengths and Ir-O-Ir bond angles 21 in addition to 
the bond-isotropic orbital overlaps observed in recent ab ini¬ 
tio results . 30 For the 7 -LIO case, the soft-spin approach with 
mixing of eigenmodes can generate the experimental phase 
using both the JKT and CZC models. However, none of the 
models can generate the experimental phase within the classi¬ 
cal approach or the soft-spin approach without linear combi¬ 
nation of eigenmodes. This may also be related to the more 


distorted nature of the 7 -LIO lattice as observed by structural 
analysis , 22 which may imply large bond-anisotropies in the 
effective spin model. These anisotropies together with the ab¬ 
sence of any symmetry on the X/Y bonds suggests that other 
symmetry-allowed, subdominant exchanges not yet accounted 
for in either the JKT or CZC model may be important in de¬ 
termining the ground state in the real material. To assess these 
possibilities, ab initio studies like DFT or quantum chemistry 
calculations may provide the necessary insights and are thus 
crucial in our understanding of 7 -LIO. 

As noted in Sec. IV A, the classical solution, as determined 
by LT or by any other method, yields a definite variational 
wavefunction and variational energy of the quantum spin- 1/2 
model. At the cost of obtaining a definite wavefunction and 
energy, the classical approach ignores all quantum effects. Al¬ 
though the classical approach may correctly generate the bro¬ 
ken symmetries deep within the long-range ordered phases, 
near phase boundaries where quantum fluctuations may be en¬ 
hanced, this approach may fail to identify the correct broken 
symmetries (if there are any) of the ground state. It may also 
be possible that the classical approach fails completely when 
quantum fluctuations are sufficiently large that the long-range 
order is destroyed. Despite these shortcomings, some quan¬ 
tum corrections to the wavefunction and energy in addition 
to the stability of the phases can be partially accounted for 
via spin-wave theory or Jastrow factors in variational Monte 
Carlo schemes. 

It was suggested 29 that the soft-spin approach may generate 
spin structures that incorporate the effects of quantum fluctu¬ 
ations. The procedure of forming linear combination of mul¬ 
tiple eigenmodes allows for construction of more complex or¬ 
dering at the same wavevector. As a trade-off, the resulting 
expectation value of the Hamiltonian using this soft-spin state 
(in both unmixed and mixed cases) cannot be interpreted as 
a variatonal energy and the obtained spin structure cannot be 
readily translated into a wavefunction. Therefore, whether the 
classical or soft-spin result is favoured in a quantum treatment 
at finite or zero temperature cannot be determined within these 
two methods alone and a comparison with a quantum treat¬ 
ment would be ideal. However, the study of such 3D frustrated 
quantum spin models is currently computationally prohibitive, 
especially when searching for incommensurate phases. There¬ 
fore, to gather more insight regarding the methods and mod¬ 
els considered thus far, we consider a ID spin chain limit of 
both the JKT and CZC models, which can be tackled by the 
DMRG technique. 


V. ID SPIN CHAIN MODEL 


The ID spin chain model related to both the JKT and CZC 
models is an extension of the spin chain model considered in 


7 



(a) DMRG - decoupled JKY (b) DMRG - decoupled CZC 



Ref. 29. It is given by 


H = ]T Ji Si ■ Sj + KS7S] + e i6 'i Sf 

(ij)eaP( 7 ) 

+ ^ Jj t ■ S r , (3) 

«ij » 


where only X and Y bonds are included in the sums, ((ij)) 
denotes second nearest neighbors (2NN), and the phase fac¬ 
tor appearing before F is defined identically to Eq. 1. When 
r/0 and J 2 = 0, we arrive at the decoupled limit of the 
JKY model: it is identical to Eq. 1 with vanishing Z bond ex¬ 
changes, resulting in decoupled chains with only X/Y bonds. 
When r = 0 and J 2 = 0, we similarly reach the decoupled 
limt of the CZC model. The 2NN Heisenberg exchange has 
been introduced to this decoupled CZC limit in Ref. 29 to 
compensate for the loss of the Z bond exchanges . 29 By study¬ 
ing these two limits with all three approaches—DMRG, clas¬ 
sical, and soft-spin analysis—we aim to compare the three 
methods and to understand the essential ingredient needed to 
generate incommensurate correlations. A more in-depth dis¬ 
cussion of the DMRG results can be found in Appendix B. 



FIG. 4. (Color online) Phase diagram of ID spin-chain model in 
the decoupled JKY limit (F < 0, J 2 = 0) and decoupled CZC 
limit (T = 0, J 2 < 0) obtained via three different methods. The 
regions enclosed by the red dotted lines indicate that LT failed. In 
a), c), and e), we studied the decoupled JKY limit with DMRG, 
classical, and soft-spin analysis respectively. Similarly, in b), d), and 
f), we simulated the decoupled CZC limit. In the decoupled JAT 
limit, all three methods predict an incommensurate phase (Sp) and an 
antiferromagnetic (AF) phase. These phases are found in analogous 
positions in parameter space compared to the 3D JKY models in 
Fig. 2. In the decoupled CZC limit, there are three commensurate 
phases: away from the stripy (ST) phase at small J 2 , a ferromagnet 
(FM) and antiferromagnet (AF) appear. The classical d) and soft- 
spin f) analyses also predict an intermediate incommensurate spiral 
(Sp) phase that does not appear in the quantum analysis. The dotted 
lines in b) are cuts for which we show the structure factor in Fig. 5; 
the structure factor peaks for 9 = ^ in a) are shown in Fig. 6. 


In the decoupled JKY limit. Fig. 4 shows that all three ap¬ 
proaches capture incommensurate correlations. This shows 
that the spiral correlations generated by the frustration be¬ 
tween the F term and the HK terms persist in the quantum 
limit. We notice that the classical and soft-spin analysis yield 
similar phase boundaries that coincide with the area where LT 
fails (enclosed by the red dotted lines), while in the quantum 
limit, the incommensurate region enlarges beyond the red dot¬ 
ted lines. On the contrary, in the decoupled CZC limit, we find 
that the quantum phase diagram obtained via DMRG does not 
contain the incommensurate correlations predicted by both the 
classical and soft-spin analysis. In particular, the commensu¬ 
rate phases are captured by all three methods (quantum, clas¬ 
sical, and soft-spin), but an intermediate spiral phase is only 
present in the classical and soft-spin analysis within the re¬ 
gion where LT fails. This suggests that the incommensurate 
correlations driven by the 2NN Heisenberg is ultimately dis¬ 
favoured by quantum effects. Taken together, these two sets 
of results suggests that finite V, and not further neighbor in¬ 
teractions, can stablize incommensurate correlations in the ID 
limit. 

One may hypothesize that the same trend will occur in the 
3D models discussed in this work: in the quantum limit, the 
Csp a phase predicted using the soft-spin approach in the CZC 
model for both lattices may be disfavoured relative to the com¬ 
mensurate phases due to the absence of T along the chain di¬ 
rection. However, since quantum effects in ID systems are 
typically greater than 3D systems, this extrapolation of results 
remains speculative. Future studies are needed to establish the 
connection between these ID decoupled chain results and the 
quantum phase diagram of the 3D models. 
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VI. SUMMARY 


We have completed the picture initated in two previous 
works by examining both the JKT and CZC models using 
both the classical and soft-spin approaches. For the hyper¬ 
honeycomb, the robustness of the JKT model in predicting 
the experimental phase regardless of the method used suggests 
that bond-isotropic interactions can explain the observed spi¬ 
ral order. On the other hand, the experimental phase of 7 -LIO 
can only be reproduced using the soft-spin analysis with mix¬ 
ing of eigenmodes using both JKT and CZC models, sug¬ 
gesting that further investigations are needed to deduce an ef¬ 
fective model from microscopic origins. In all cases where the 
experimental phase was predicted, a dominant ferromagnetic 
Kitaev exchange exists as a common thread. We conclude 
by examining the decoupled limit of both the JKT and CZC 
models and observed that a finite 1 was needed to ensure that 


the incommensurate correlations predicted in both the classi¬ 
cal and soft-spin methods persist in the quantum limit. Further 
exploration extending these findings to higher dimensions are 
warranted in order to assess the relevance of these results to 
the 3D honeycomb iridates. 
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Appendix A: Magnetic basis vectors 

The irreducible representations and magenetic basis vectors 
for a magnetic structure with ordering wavevector q = (h00 ) 
in the orthorhombic unit cell of the hyper-honeycomb lattice 
are given by 32 


Irreducible 

Basis 

Representation 

Vectors 

IT 

Fa, G b , A c 

r 2 

C a , A b , G c 

r 3 

Ga, F b , C c 

r 4 

A*, Cb, Fa 


In the sublattice basis given in Fig. la, the magnetic basis 
vectors are given by 32 


F = 


’ 1 " 


’ 1 " 


’ 1 ’ 


’ 1 " 

1 


1 


-1 


-1 


,c = 


,A = 


,G = 


1 


-1 


-1 


1 

1 


-1 


1 


-1 


(Al) 


With a stmcture specified by the basis vectors v = 
( iv a , iVb, v c ), the moments at each site is given by ’ 2 

S n (r) = aS a v a (n) sin q ■ r + hSbVb (n) sin q ■ r 

+ cS c v c (n) cos q- r, (A2) 


where a, b , c are the unit vectors of the orthorhombic lattice 
vectors, S n (r ) is the moment at position r which has sublat¬ 
tice index n, and .S', are the magnitudes of each component of 
the moment. Since all phases predicted in this work have rel¬ 
ative phases between components given by v = ( iv a , iVb, v c ), 
we adopt a shorthand notation by dropping the relative imag¬ 
inary unit i: v = (v a , Vb, v c ). 

For the stripy-honeycomb, there are 8 sublattices in the unit 
cell as opposed to 4 for the hyper-honeycomb. We break up 
those 8 sublattices into two sets of 4—the unprimed 1 — 4 sub¬ 
lattices and the primed sublattices 1' — 4' as seen in Fig. lb. 33 
The irreducible representations and magnetic basis vectors for 
wavevector hOO are given by 33 




Ji 


Ji 


(a) Ji = 0.33 


(b) Ji = 0.25 


FIG. 5. Decoupled CZC model: specific cuts of the static structure 
factor I(q) = for q = 0 (FM), tt/ 2 (ST) and n (AF) 

are shown in (a-b) for N = 128 (solid lines) and N = 96 (dashed 
lines). Along with these we shown the maximum of I(q), which 
always belongs to one of these three choices. 



FIG. 6 . Decoupled JKY model: location of the maximum static 
structure factor (i.e. ma x q [I(q) = I MM (9)]) for 6 = 37 t/4 as a 

function of <j> for the N = 128 case. The antiferromagnetic phase 
for (j) > 77 t /4 has the maximum peak located at ±tt. Below this 
value of (j>, incommensurate correlations develop and persist down to 

< f > = 3tt/2- 


Appendix B: ID spin chain results 


Irreducible 

Basis 

Representation 

Vectors 

IT 

C a , A„, Ga 

r 2 

F a , G b , A c 

r 3 

A a , C b , F c 

r 4 

G a , F b , C c 


The magnetic basis vectors are also specified by Eq. Al. 
The structure given by 

v = i(v a ,-v' a ),i(v b , -v' b ), (v c , v' c ) (A3) 

is also specified by Eq. A2 but with a relative sign on the a 
and b components between the primed and unprimed sublat¬ 
tices. Since all magnetic structures predicted in this work has 
the relative phases given by Eq. A3, we adopt the shorthand 
notation where v = ( v a , Vb, v c ). 


To determine the ground state phase diagram of Eq. (3), we 
used the DMRG as implemented in the ITensor framework. 
We consider a chain of N < 128 sites with open boundary 
conditions. In the decoupled JAT limit where J 2 = 0, we 
parameterize our phase diagram in the same manner as per Eq. 
2. In the decoupled CZC limit where T = 0, we fix the energy 
scale K = — 1 such that there is a two-dimensional parameter 
space (Ji, J 2 ) to explore. For each point in parameter space 
we perform 10-20 sweeps, keeping up to 1000 states resulting 
in a maximum truncation error of ~ 10 9 . For small system 
sizes (N < 22) these results were checked for consistency 
against exact diagonalization. 

To characterize the phases we consider two diagnostics: 
derivatives of the ground state energy and the static structure 
factor. Sharp features in the energy derivatives signal phase 
boundaries while peaks in the diagonal components of the 
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static structure factor 

Iw( q )=±Y, eiq ' {i ~ j) (S?S»), (Bl) 

ij 

signal the ordering wave-vector. 

For the decoupled CZC model, we find that throughout the 
entire phase diagram the maxima of the structure factor I(q) 
appear only at wave-vectors q = 0, 7 r /2 or tt corresponding 
to a ferromagnet (FM -xy), antiferromagnet (AF-z) or stripy 
(ST) phase. The isotropic structure factor I = I^(q) for 

each of these wave-vectors is shown in Fig. 5 for both N = 96 
and N = 128. The ferromagnet and stripy phases primarily 
have correlations in the x and y directions with the maxima in 
appearing in I xx or I vv , while the anti-ferromagnet is oriented 
in the z direction, with the corresponding maximum in I zz . 

The gross features of this phase diagram can be understood 
easily. First we have two exactly solvable points present in 
the J 2 = 0 limit; the Kitaev chain at .J t = 0 and a point dual 
to a Heisenberg ferromagnet at J\ = 1/2. As in the two- 
dimensional case, a stripy phase extends between these two 
solvable points. This is stable to finite J 2 , but ultimately gives 
way to either the VM-xy or the AF-z phases. These can be un¬ 
derstood by considering the large J 2 limit where the even and 
odd sites decouple into two ferromagnetic chains. The relative 
orientation of the two sublattices is then determined by the I\ 
and Ji interactions, with K < 0 favouring ferromagnetism 
and J\ > 0 favouring anti-ferromagnetism. 


The anti-ferromagnetic phase for J\ > 0.3 can also be 
understood by considering the solvable point at (Ji, J 2 ) = 
(1/2,0). This point is dual to the ferromagnetic Heisen¬ 
berg point via a four-sublattice transformation, and hence pos¬ 
sesses SU(2) spin rotation symmetry in the rotated basis. One 
finds this continuous degeneracy is broken by finite J 2 , se¬ 
lecting the AF-z state as the ground state, which is a simple 
product state. Since this product state is the ground state of the 
ferromagnetic J 2 interaction and the (Ji, K) = (1/2, —1) in¬ 
teraction independently, it follows that this product state is the 
exact ground state of the (J-|, K, ,/ 2 ) = (1/2, —1, J 2 ) line. As 
Ji decreases, the AF correlations are preserved but the ground 
state is no longer a simple product state. 

For the decoupled JKT model, we find one boundary 
across the phase diagram as indicated by energy derivatives. 
The maxima of the structure factor of the phase connected to 
the antiferromagnetic Heisenberg point is located at ±t r, in¬ 
dicating AF correlations. As we move across the boundary 
to the adjacent region, peaks in the structure factor develop at 
incommensurate wavevectors while the peak at ± 7 r vanishes. 
This incommensurate peak’s position varies as we approach 
the Ji = 0 limit. The stripy phase located in the 1 ’ = 0 
limit does not extend into the finite F region, in contrast to 
the effects of J 2 in the decoupled CZC model. We illustrate 
the position of the maximum of the static structure factor for 
6 = 37 t/4 as a function of (j> in Fig. 6. 


